Shot noise from action correlations 
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We consider universal shot noise in ballistic chaotic cavities from a semiclassical point of view 
and show that it is due to action correlations within certain groups of classical trajectories. Using 
quantum graphs as a model system we sum these trajectories analytically and find agreement with 
random-matrix theory. Unlike all action correlations which have been considered before, the corre- 
lations relevant for shot noise involve four trajectories and do not depend on the presence of any 
symmetry. 
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One of the most prominent methods to describe spec- 
tral and transport properties of ballistic quantum sys- 
tems with classically chaotic dynamics relies on the semi- 
classical representation of the Green's function in terms 
of classical trajectories. For closed systems this leads to 
Gutzwiller's trace formula [l|, which expresses the oscil- 
lating density of states as sum over periodic orbits. For 
open systems the semiclassical theory of chaotic scatter- 
ing and in particular its applications to electronic 
transport through mesoscopic devices 4], are based on 
this approach. The power of this method as compared, 
e. g., to random-matrix theory (RMT) 0, @ is its po- 
tential to account for system-specific details. Its main 
drawback lies in the difficulty to handle the resulting 
sums over huge sets of classical trajectories. Until re- 
cently the only method to deal with this problem was 
Berry's diagonal approximation 0, which neglects any 
nontrivial correlations between the trajectories. As a 
consequence many interesting phenomena such as weak 
localization, universality in spectral statistics and in con- 
ductance fluctuations, or the supression of shot noise can- 
not be described properly. Although the role of correla- 
tions between the actions of classical orbits has been ap- 
preciated for a long time [8-10], there is only one special 
case where they can be accounted for explicitly: Sieber 
and Richter recently calculated the leading order weak lo- 
calization corrections from correlations between specific 
orbit pairs II 1| . This result stimulated further intense 
research |l2l 1 1 3| and it is clearly a very promising ap- 
proach. However, as mentioned above, weak localization 
is just one besides a variety of other phenomena for which 
action correlations are relevant as well, but cannot be ac- 
counted for by the orbit pairs considered in [11-13]. 

In this paper we address shot noise in ballistic meso- 
scopic conductors, which is an important source of (ex- 
perimentally accessible) information about the dynamics 
in such systems [18-21]. We identify the action corre- 
lations which are responsible for shot noise and explain 
how due to these correlations the universal result of RMT 
[5, 14-20] can be recovered from a single system without 
ensemble average. The relevant correlations are funda- 



mentally different from all those considered previously 
[8-13] because they involve four instead of just two clas- 
sical trajectories. Moreover they do not depend on the 
presence of symmetries. We emphasize that action cor- 
relations are no small correction. They are needed to 
understand shot noise to leading order. 
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FIG. 1: (a) A chaotic cavity with two attached waveguides 
and a classical trajectory contributing to conductance and 
shot noise, (b) The quantum graph used to model this situa- 
tion. The transversal modes of the waveguides correspond to 
the infinite leads attached to the graph (bold lines), while the 
internal system is represented by many finite bonds. 

We will perform all explicit calculations in a specific 
model system, the quantum graph of Fig. \T]p. Quan- 
tum graphs (networks) have a long record as models for 
electronic transport (see (2^| and Refs. therein). Since 
the pioneering work of Kottos and Smilansky they are 
also established in quantum chaos [23-25, 12]. They are 
particularly suitable for our purpose as the representa- 
tion in terms of classical trajectories is exact and also 
the analogue of action correlations amounts to exact de- 
generacies. Previous studies showed that despite these 
analytical simplifications the mechanism and the role of 
correlations between classical trajectories are equivalent 
to other systems such as billiards |l^ . 

Consider a chaotic cavity with two attached wave- 
guides supporting N±, N2 transversal modes, respectively 
(Fig. [TJi). For small bias voltage and temperature, negli- 
gible electron interactions and fully coherent dynamics all 
information about electron transport through this system 
is contained in the scattering matrix at the Fermi energy 
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Shot noise represents temporal current fluctuations due 
to the discreteness of the electron charge [26j. At zero 
temperature the average power of the noise can be ex- 
pressed in terms of the transmission matrix t as 0, |27} 

P = 2e\V\G Trttt(l-ftt), (2) 

while the assumption of uncorrelated electrons yields 
P Poisson = 2e\V\G. Here G = G Trtt f denotes the 
conductance, e is the electron charge, V the voltage and 
Go = 2e 2 /h the conductance quantum. RMT yields Q 

(P)rmt = 2e\V\ Go NlN*/{N x + N 2 ) 3 (3) 

and for the conductance (G)rmt = Go A^iiVa/(iVi 
Each result comes with a weak localization correction 
which is small {(SP) < (P), (5G) < (G)) for large N u 
N2 and will not be considered here. For N% = N2 the 
Fano factor F — (P)/(Pp isson) = 1/4 is obtained. 

We will show that Eq. ® can be recovered semiclas- 
sically. To this end, the element t n2ni of the trans- 
mission matrix is expressed as a sum over all classical 
trajectories connecting the incoming mode 1 < ni < 
Ni with the outgoing mode 1 < n 2 < N2 [3 

t n2ni =Y,A P e iS */ n , (4) 
p 

where S p denotes the classical action and A p is an ampli- 
tude related to the stability of the trajectory. While for 
the conductance we have to evaluate a sum over pairs of 
trajectories 

Ni N 2 

<Trttt )= £ £ £ ^<eW-*>>, (5) 

ni—l n,2 — l pq:n\ — m2 

the shot noise involves also a term combining four clas- 
sical paths 

(Tr^t)^ = £ £ ]Tvi^: 

mini m2ri2 pgrs 

x / e i/ft(Sp-S 9 +Sr-^)^ (g) 

Here, the trajectories p, q, r, s connect two incoming to 
two outgoing modes as shown in Fig. [2t. As Eqs. JHJ), © 
describe one particular system rather than an ensemble, 
the average (•) is to be taken over an energy window. 
It should be small enough to keep the classical dynam- 
ics and in particular the amplitudes A p essentially un- 
changed. Nevertheless, in the semiclassical limit Ti — > 
the phase factor is rapidly oscillating and only those orbit 
combinations survive the averaging for which the action 
changes are correlated. 

In particular, setting p = q in Eq. 101 the phase drops 
out and we are left with a sum over classical probabili- 
ties |^4 P | 2 - This is the diagonal approximation 0. Pro- 
vided that the dwell time inside the open cavity is large 



compared to the time needed for equidistribution over 
the available phase space, the probability is the same 
for all outgoing modes and (G)rmt is exactly recovered 
0. Hence, the contribution from other pairs of corre- 
lated trajectories that might exist must vanish although 
no explicit demonstration of this fact has been given up 
to now. In presence of time-reversal symmetry the above 
remains valid to leading order in the mode number N. 




FIG. 2: (a) shows schematically four classical trajectories 
p, q, r, s connecting two incoming modes mi, ni to two outgo- 
ing modes mi,ni such that the diagram contributes to the 
semiclassical approximation of Eq. (J2J. A contribution to 
Eq. (JSJ results only if the trajectories are correlated such that 
the action difference between p, r (full lines) and q, s (dashed 
lines) remains small for variying energy, (b) and (c) shows 
the simplest configurations where this is the case: two trajec- 
tories are pairwise equal (diagonal approximation) . (d) shows 
the configuration which completely accounts for the universal 
shot noise in leading order. 

For shot noise the diagonal approximation has two ana- 
logues (Fig. 03,0): we can have (i) p = q, r = s for 
mi = rii and (ii) p = s, r = q for m,2 = ri2- In both cases 
no phases are left in Eq. © and the remaining summa- 
tion is over two independent classical trajectories p, r 
with the only constraint that they begin or end at the 
same mode, respectively. Proceeding as in the diagonal 
approximation to the conductance we obtain [3l| 

m i—l ni2 ,n-2 — l pr 

rn 1 , n 1 — 1 mo — 1 pr K ' 

Combining these two results we find 

(Tr (tit) 2 ) diag = (Trttt) diag = N^/i^ + N 2 ) (9) 

and according to Eq. this means that within the di- 
agonal approximation there is no shot noise. This is no 
surprise: the diagonal approximation reduces the quan- 
tum to the classical problem and since classical dynam- 
ics is deterministic there is no uncertainty if an incoming 
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electron is transmitted or not and hence no noise |l8tl29| . 
On the other hand Eq. © is quite remarkable as it means 
that within the semiclassical approximation shot noise is 
entirely due to nontrivial correlations between different 
trajectories. 

What is the general mechanism for such correlations? 
Previous research [9-13] showed that pairs of trajecto- 
ries have correlated actions if they explore the same (or 
symmetry-related) parts of phase space with a different 
itinerary. In terms of symbolic dynamics the code words 
for the two orbits are composed of the same sequences, 
in permuted order. The analogy to dia gram matic per- 
turbation theory and some recent results [llj, |l2j suggest 
further that the importance of the correlations decreases 
with growing number of sequences needed to represent 
the code of the trajectories: in the diagonal approxima- 
tion to the conductance the codes are equal and the result 
is correct to leading order, orbit pairs composed of two 
loops give the next-to-leading order correction etc. 

In the case of shot noise we have exhausted the di- 
agonal approximation and consider therefore trajectories 
p, q, r, s whose codes can all be represented in terms of 
two subsequences. Inspection shows that the only option 
is the diagram of Fig. 01. Indeed the phase in Eq. © 
will almost vanish for such contributions since the combi- 
nation of p, r (full lines) almost coincides with the com- 
bination of q, s (dashed lines) such that the respective 
actions cancel. A remaining small total action difference 
comes from the different behaviour of the trajectories in- 
side the crossing region. In this respect the correlated 
trajectories of Fig. EJl are very similar to those giving 
rise to weak localization effects [11-13], i. e. the methods 
developed there for various specific systems should allow 
for a straightforward generalizion to shot noise. 

In the remainder of this paper we treat one of those 
systems explicitly, namely the quantum graph shown in 
Fig. 03. The closed version of this graph consists of a 
central vertex with valency B and b = 1 . . . B attached 
bonds with incommensurate lengths L b . Following the 
standard quantization j23l | the dynamics of a particle 
with wavenumber k = (2mE) 1 ^ 2 /h is represented by a 
B x B bond-scattering matrix Yj\,b'{k) = crbb'& 2lkLb con- 
taining energy-dependent phases 2kL b from the free mo- 
tion on the bonds and complex amplitudes describing the 
scattering at the central vertex. A basic requirement is 
unitarity (current conservation) 

B 

<W = 2J <7 PaO-0 a , , (10) 

13=1 

otherwise a can be chosen according to the physical sit- 
uation. For simplicity we set 

a w = e W/s /V5 (11) 

such that all classical transition probabilities are equal 
|<7bb'| 2 = 1/B. This model was first considered by Tanner 



[25j who showed numerically that its specral statistics fol- 
lows RMT. We open the graph by extending N — Ni+N 2 
bonds to infinity and model a two-channel geometry by 
considering N± (N 2 ) of these leads as the modes in the 
left (right) contact (Fig. Hb). For fixed N%, N 2 we will 
consider the limit B — ► 00 in order to meet the condi- 
tion of a long dwell time which was already mentioned 
in connection with the diagonal approximation. Within 
the leading-order in B we also neglect lower-order correc- 
tions in the mode numbers N\, N 2 in order to compare 
our result to Eq. J3J. 

For a graph the N x N unitary scattering matrix 
Eq. (0) can be expressed in terms of subblocks of the 
bond-scattering matrix £ via S = ^cc + ^cg{I ~ 
Y,gg)~ 1 T,g C [23] . C = C\ U C 2 denotes here the set of 
N = Ni + N 2 leads and Q comprises the B — N bonds 
inside the graph. Expanding the Greens function of the 
internal part (1 — Yigg)~ l into a geometric series we arrive 
at Eq. Q which is in the case of graphs an identity rather 
than a semiclassical approximation. The sum is over all 
trajectories (=bond sequences) p = [n\p\ . . . ptn 2 ] con- 
necting the lead n\ £ C± to the lead n 2 S C 2 via an 
arbitrary number t > of internal bonds pj € Q . The 
action is related to the total length of the trajectory, 
Sp/h = kL pi where L p = X^j=i L Pr Finally the ampli- 
tude is given as A p = cr niPl a PlP2 . . . <r Pt n 2 such that the 
classical probability of the trajectory is \A p \ 2 = 1/B t+1 . 
Consequently, we have 

p t=0 

from which we do indeed recover the diagonal approxi- 
mation Eqs. 0-0- 

Next we consider trajectories p, q, r, s which are com- 
posed of four sequences a, 6, c, d as shown in Fig. 01 and 
assume that each of these sequences has a length t > 1. 
In order to avoid overcounting we have to ensure that for 
any given set p, q, r, s the definition of a, 6, c, d is unique. 
Potential problems arise if all four trajectories coincide 
in the crossing region for one (or more) steps: p = [(17c], 
q = [67c], r = [tryd\, s = [a7d]. It is a matter of taste if 
7 in this situations is considered as part of a and b or of 
c and d. We use the first representation and enforce it 
by the restriction a ^ d\ (The subscripts i/f are used for 
the initial/final bond in a, 6, c, d.). 

Returning to Eq. JHJ we note that the actions of p, r 
and o, s cancel exactly such that the phase factor is ab- 
sent |32(. As in Eq. I|12|l we can perform the summation 
over all internal bonds of the subsequences a,b,c,d and 
also over the leads mi, n\, m 2l n 2 . Only the amplitudes 
from transitions right at the intersection of a, b, c, d do 
not combine into classical probabilities and must be con- 
sidered explicitly. We obtain 

(4) N?N$ \ - * » n „, 

l abcd ~ N 4 1^ a a 1Cl <J b{Ci cr btdi a aidi (13) 

G[f6f ,Ci^dj^Q 
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E ~ E E j^m^tic^hdi^di 
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(15) 



To perform this calculation we have repeatedly used 
Eq. HI L)| l in a form which allows to transfer summations 
from the graph Q to the leads L 
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Moreover we used that according to Eq. I|ll[) we have 
\ (J a{c i o'b t c i a b{d i o'atd i I = & which means that the two 
sums in the second line of Eq. I|14(l yield only a negli- 
gible correction 0(B~ 1 ) because the number of terms 
are N 3 (N - 1) and N(N — 1)(B — N), respectively. It 
is the first line in Eq. I|14[l which gives the dominant 
contribution. With exactly the same methods we can 
consider the contribution from special cases of the di- 
agram in Fig. [2i where the length of one of the sub- 
sequences a, b, c, d vanishes. We find that for vanish- 
ing a or b we get = t^ d = -t^L while for van- 
ishing c or d no contribution results. Finally we have 



(Tr(ttt) 5 



t {i) 

abed 



b bcd 



(4) 
acd 



-NfN$/N 3 and sub- 



stitution of this result into Eq. shows that we have 
indeed reproduced the RMT result from correlated clas- 
sical trajectories. 

We have checked that our result remains unchanged if 
we substitute in Eq. © Tr it* (1- tit) = TrttW. Quan- 
tum mechanically this is just a consequence of unitarity, 
but within the semiclassical approach it is a nontrivial 
result since entirely different trajectories contribute and 
unitarity is restored only if all relevant correlations be- 
tween them are properly accounted for. 

We are grateful to O. Agam and H. Schomerus for 
stimulating our interest in the problem. 
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